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I. INTRODUCTION 

Newly discovered iron-pnictide superconductors have 
a number of uncommon properties. The subject of this 
paper are two such properties: (a) The specific heat jump 
AC is proportional to T 3 as demonstrated on "122" se- 
ries of Ba(Fei_ 2: 00^)2 As2 and Ba(Fei_ :c Ni x )2As2. 1 This 
behavior, according to Ref. 2, cannot be explained within 
the "realm of conventional BCS theory" . Similar behav- 
ior is recorded in 122 crystals with Ba substituted par- 
tially with K and with Fe substituted with Pd, Rh, 3 and 
Co-Cu. (b) Slopes of the upper critical field dH c2 /dT at 
T c are proportional to T c across both 1111 and 122 series. 

It is shown below that both scalings can be under- 
stood within the weak-coupling BCS model provided a 
strong pair breaking is present in these materials. In fact, 
these features should also be present in conventional su- 
perconductors with magnetic impurities as discussed by 
Abrikosov and Gor'kov (AG) in their seminal work on 
the pair breaking for the nearly critical concentration of 
these impurities when T c <C T c q, the critical temperature 
of clean material. 4 AG had considered isotropic materi- 
als with a spherical Fermi surface and the s-wave order 
parameter constant along this surface. The symmetry of 
the order parameter in multi-band pnictides is not yet 
determined with certainty; however, many favor the ±s 
structure. 5 ' 6 The critical temperature in materials with 
a strongly anisotropic order parameter is suppressed not 
only by scattering breaking the time reversal symmetry 
(e.g., the spin-flip); in fact, any scattering reduces T c . 
The term "pair-breaking" is used here in a broad sense 
for any process suppressing T c . It is shown below that 
both features, dH c2 /dT oc T c and AC oc T c 3 , follow from 
the assumption that the "pair-breaking in a broad sense" 
is strong. 

Below, the linearized Ginzburg-Landau (GL) equation 
and the energy near T c are derived within the weak cou- 
pling scheme (that allows one to evaluate dH C 2/dT and 
AC at T c ) for an arbitrary anisotropy of the order pa- 
rameter A and of the Fermi surface in the presence of 
pair-breaking. The text is focussed on the situation when 
the average (A) over the F-surface is close to zero that 
presumably is the case of pnictides. 5 ' 6 Comparison with 
the data available concludes the paper. 

Perhaps, the simplest for our purpose is the Eilen- 
berger quasiclassical formulation of the Gor'kov's theory 



that holds for a general anisotropic F-surface and for any 
gap symmetry: 7 

vUf = 2Ag -2u;f+ — (f) - -L(g) , (1) 

T_ 7 + 

-vU*f+ = 2A*g - 2ojf+ + -L{f+) - f —{g) , (2) 



g 2 = 



(3) 



A(r,k F ) = 2nTN(0)Y / (v(k F ,k' F )f(v , ,r,uj)) ; (4) 



ui>0 

Here, v is the Fermi velocity, IT = V + 2iriA/4>o, (f> Q is 
the flux quantum. A(r, hp) is the order parameter that 
in general depends on the position k F at the F-surface of 
other than the isotropic s-wave symmetry. The functions 
f(r,v,uj), /+, and g originate from Gor'kov's Green's 
functions integrated over the energy variable near the F- 
surface. Further, N(0) is the total density of states at 
the Fermi level per one spin; the Matsubara frequencies 
are u> — ttT(2ti + 1) with an integer n and h = feg = 1. 
The averages over the F-surface are shown as (...). 

The scattering in the Born approximation is charac- 
terized by two scattering times, the transport scattering 
time r responsible for conductivity in the normal state, 
and r m for spin-flip processes: 
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— = - ± — 

T± T T m 



(5) 



The strong scattering in unitary limit is not considered 
here. Commonly, the scattering is characterized by two 
parameters 
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P : 



2nT r T 



and 



2ttT c t„ 



(6) 



or equivalently by p± — p ± p m . This is of course a sim- 
plification; for multi-band F-surfaces one may need more 
parameters for various intra- and inter-band processes. 
This and other simplifying assumptions notwithstanding, 
the model employed is amenable for analytic work and 
may prove useful. 

Long experience in dealing with pair-breaking effects 
has shown that the formal AG scheme in fact applies to 
various situations with different causes for the pair break- 
ing, not necessarily the AG spin-flip scattering. 8 In each 
particular situation, the parameter p m must be properly 



2 



defined. Without specifying the pair breaking mecha- 
nism in materials of interest here we apply below the AG 
approach to show that the pair breaking accounts for 
experimental data on slopes of H c2 at T c and for quite 
unusual dependence of the specific heat jump on T c . 

Commonly, the effective coupling V is assumed factor- 
izable, V(k F ,k F ') = V Q(k F ) Q(k F ' ). 9 One then looks 
for the order parameter in the form: 



A(r, T; k F ) = *(r, T) VL{k F ) 



(7) 



Our notation is motivated by the fact that so defined 
^(r,T) enters the Ginzburg-Landau (GL) theory near 
T c . The function £l(k F ), which describes the variation of 
A along the F-surface, is conveniently normalized: 



1. 



(8) 



Then, the self-consistency equation (4) takes the form: 



(r, T) = 2TTTN(0)V a £ (n(k F )f(k F , r, w) 



UJ>0 



(9) 



The assumption of a factorizable potential is quite re- 
strictive as far as complicated F-surfaces and interac- 
tions are concerned. E.g., within a two-band scheme with 
four independent coupling constants , the factorizable 
model implies V11V22 — V12V21 = 0. 

Instead of dealing with the effective microscopic 
electron-electron interaction V and with the energy scale 
u>d, one can use within the weak coupling scheme the 
critical temperature T c q (of the hypothetic clean mate- 
rial) utilizing the identity 



^(0)^0 



In- 



l c0 



t—^ in 



(10) 



UJ>0 



which is equivalent to the BCS relation Ao(0) = 
7rT c0 e~ 7 = 2uj d exp(— 1/N(0)V ); 7 is the Euler constant. 
Substitute Eq. (10) in Eq. (9) and replace ujd with infin- 
ity due to fast convergence: 



LU>0 



0/ 



(11) 



II. GL DOMAIN AND T c (r,r m ) 



Near T c , g = 1 - //+/2 and Eq. (1) reads: 
1 (f) 

2 \ 2t_) 2t+ 



Here. 



lo+ = u) + 1/2t+ , 



(12) 



(13) 



and the terms on the RHS are arranged according to the 
their order in powers of St = 1 — T/T c : the terms on the 
upper line are of the order St 1 / 2 whereas on the lower 
line - St 3 / 2 . Note that on the LHS, Uf ~ //£ ~ St. 

We look for the solution f = f\ + fi + ... where /1 ~ 
St 1 ' 2 and f 2 - St. Hence, we have in the lowest order: 



= A - oj+h 



(A) 

2r_ 



(14) 



Taking the average over the Fermi surface we obtain 

= (A)/u m , uj m =uj + l/r m (15) 
(note the difference in definitions of lo + and u> m ). Hence: 



h = 4 ( a 



(A) 
2r-Uj ri 

Comparing terms of the order St, we get 

!/,) = -«=», 



and 



h = - 



1 



2w + w„ 



2uj r , 



-vU A 



(A) 
2r_w K 



(16) 



(17) 



(18) 



Evaluation of higher order corrections for arbitrary A 
anisotropy is increasingly cumbersome unlike the case 
(A) = for which one finds: 



2u>j r ui m 



A 



(A 2 ) 
2t+uo^ 



(19) 



The critical temperature of materials with anisotropic 
order parameter is suppressed by scattering. In zero 
field, all quantities are coordinate independent; besides, 
as T — > T c , g — ► 1. Therefore, we can utilize / of Eq. (16) 
in the lowest order to obtain for T r : 



2irT c T c 



OO 

w>0 



1 



2w" l Wc 



-), (20) 



where the subscript c is to denote that w's are taken at 
T c . This generalization of the well-known AG result gives 
the T c suppression for any (Born) scattering for arbitrary 
symmetry of the order parameter; it has originally been 
obtained by Openov. 10 The sums here are expressed in 
terms of di-gamma functions: 



In 



T, 



c0 



1 



- w 



1> 



1> 



Pn 



(21) 



If T c — > 0, one can use asymptotic expansion ip(z) = 
In z — l/2z for large arguments since p, p m — ► 00. The 
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leading term then gives that T c = when scattering times 
satisfy the relation: 



1 / r m y-<n) 2 _ A (0) 
r m \2t+J 2 



(22) 



Here, A o (0) = nT c0 e~ 7 is the zero temperature gap of 
the (hypothetic) scattering free material. Clearly, this 
reduces to the AG critical rate l/r TO = A o (0)/2 for 
isotropic order parameters. If (f2) = (as, e.g., for the d- 
wave) , we have the critical combined rate: l/r+ = A (0) . 

In the absence of spin- flip scattering (r m — > oo) the 
LHS is zero and Eq. (22) has no solutions for r, i.e., T c 
does not turn zero for any r. However, a finite r at which 
T c = does exists for any finite r m . One can show that 
near the critical value rf^ 1 = A (0)(T m /2)< n > 2 , the 
critical temperature behaves similarly to the AG gapless 
case, T c ^(r + -r+ lt f/\ 

Combining Eqs. (11) and (20) one excludes the unphys- 
ical T c0 : 



U1>0 



2ttT 



tut 



+ 



*(ft) 2 



nf)\ (23) 



where t = T/T c . 



Eqs. (16), (18), and (19) for / and the self-consistency 
equation. For the case of exclusively transport scattering 
(r m = oo), the GL equations have been derived in Ref. 13. 
It is done below taking a finite r m . 

To write the self consistency Eq. (25) near T c one has 
to express (O/) with the help of Eq. (12). To this end, 

one applies (ft/(j + ...^ to (12) keeping terms up to the 
order St: 



0J+ 



2t_cj + \2u> + 



-vUf 



(26) 



and substitutes the result to Eq. (25): 



= E 

UJ>0 



*(ft) 5 



/ n 



2iw mc tj +c T_ 2o)+r_ 

(27) 

Since we are expanding in powers of \[S~i, the distinction 
between, e.g., u> c and u> = u> c (l — St) is relevant. 

When substituting here / = /i + fa of Eqs. (16) and 
(18) note that (fi vIIA) = because the angular depen- 
dence of fl (the symmetry of A) has nothing to do with 
that of the vector IIA. We then obtain: 



III. THE CASE T c <C T c0 

Situations of interested here are of T c strongly sup- 
pressed relative to T c0 (similar to the gapless supercon- 
ductivity of AG, but not necessarily the same). It is con- 
venient for this purpose to rearrange Eq. (23) by adding 
and subtracting ^ /lo + under the sum. We transform: 



OO 

= E 



n=>0 



-lnf- 



n + l/2 + p+/2 

£l + I 

2t 2, 
l-t 2 



n+l/2 + p+/2t 



2 2 



6pl 



(24) 



The parameter p + is large if T c — > and one can use large 
arguments asymptotics of the di-gamma functions. Com- 
bining Eqs. (23) and (24) we obtain the self-consistency 
equation in the form: 



^{l-t 2 ) 
\2-nTp\ 



w>0 



*(0) 2 



- W 



(25) 



V5t 
6nTp 2 + 



^ V2t^I 



w>0 

1 

4^7 



(«n) 2 * (^ 2 + 



- i 



i 



2t_w„ 



2aj rm T- 



(28) 



Note that the LHS and the term at the lower line of this 
equation are of the order St 3 / 2 ; for this reason all lj's in 
this term are taken at T c . Besides, the round brackets at 
the upper line of the RHS are easily shown to turn zero 
at t = 1. Expanding the bracketed expression in powers 
of St and keeping only the first term one obtains: 



A^St = -B ik UiU k ^ 



(29) 



with 
A 



QttT cP % 

(O) 2 ^ [LO l _ 



2t_ 



(2uj m + uj + )lo 



w>0 



2T-U 2 m u\ 



Bj 



4 ^ w, 

UJ>0 + 



2uo m T- 



(30) 
(31) 



where all o/s are at T c and the subscript c is omitted. 
This is, in fact, the linearized anisotropic GL equation 



A. Linearized GL equation and the coherence 
length 

The GL equations arc obtained by utilizing smallncss 
of A/ iv and of vHA/iu 2 near T c . Hence, one can use 



-(£ 2 ) ife iLn fe * = *. 

with anisotropic coherence length given by 
(Z 2 )ik = B ik /A8t. 



(32) 



(33) 
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All sums in Eqs. (30) and (36) are expressed in terms of 
poly-gamma functions of large parameters p± . Keeping 
the leading terms we obtain: 



1 



Bik — 



+ 



2ttT c 

(n)(nviv k )Tl 



7TT C(0 _ 



hn- P+ 



2p n 



In 



P+ P-{ 2 P+ + P-) 



2p r , 



2p\ 



(34) 



(35) 



B. Materials with (fi) = near T c 

This corresponds, e.g., to the d-wave symmetry. 
Within a two-band model for iron-pnictides the order 
parameter has a ±s structure, so that (A) <C |A„ laa; |. 6 
One then expects the model with (O) = to hold at least 
qualitatively for iron-pnictides. 

If (O) = 0, A and B are simplified: 



1 



6ttT, 



(Vt 2 V t V k )Tl 



We then have: 



cP+ 



2ttT c 



A-K 2 T 2 5t 



(36) 



(37) 



For the d-wave order parameter and isotropic 2D Fermi 
surface, O = v / 2cos2<p and (Q 2 v 2 ) = v 2 /2: 



2„,2 



e = 



3h 2 v 
8ir 2 T 2 5t 



(38) 



This result has been obtained in Ref. 11 for a clean d- 
wave with a strongly suppressed T c . 

For a uniaxial material, the slope of the upper critical 
field along the c direction near T c is given by 



dH, 



c2,c 



dT 



2n<j) k 2 B 

■ih 2 {n 2 v 2 a ) 



(39) 



(in common units). Although the scattering and pair- 
breaking parameters do not enter this result explicitly, 
they affect i? C 2,c and its slope via T c (p + ). One readily 
obtains for the other principal direction: 



dH, 



c2,ab 



27TC/f>0 fc| 



<it 3h 2 ^/(n 2 v 2 )(n 2 v 2 ) 



(40) 



It is worth recalling that in isotropic materials with 
the standard s-wave order parameter the slope H' c2 oc T c 
in the clean limit (because H C 2 oc l/£ 2 oc T 2 ) whereas for 
the dirty case H' c2 is T c independent (H C 2 oc l/££ oc T c , I 
is the mean- free path). The propotionality H' c2 to T c is a 
property of the AG gapless state. In our case, the result 
(39) is obtained for a strong pair-breaking in materials 
with anisotropic order parameter. 



Note also that even without magnetic scatterers, in 
materials with = and p + 3> 1 , the superconductiv- 
ity becomes "gapless" in a sense that the total density of 
states at the Fermi level is not zero. As in the AG case, 
if T c — > 0, the superconductivity is weak at all tempera- 
tures, i.e., / << 1 and g = l-.f 2 /2 = l-A 2 /2w+2 hi the 
whole domain < T < T c . Then the energy dependence 
of the total density of states N(e) — N(0) Beg(hcj — > ie) 
reads: 



N(e) 



2A 2 - 2 



1 



N(o) + (TTW V = 2T+e - (41) 

Hence, at zero energy, N(e) has a non-zero minimum, 
whereas the maximum of N(e) is reached at e rn = 
\/3/2t+ (not at A). Therefore, the ratio of the "ap- 
parent gap" e m to T c should vary as l/T c . Since only 
the total density of states is non-zero, this does not ex- 
clude possibility to have gapped and gapless patches on 
the F-surface. 



IV. THE SPECIFIC HEAT JUMP 

Eilcnberger equations (1) and (11) in zero field can be 
obtained minimizing the functional 7 



T = N(0) 



4^ In- 



T 



2.TW 



5<2 

tkij 



2A/ + 2w( 5 -l) 



LU>0 

/</) , (9(9) -1) 



2t~ 



2t+ 



(42) 



(43) 



The function g here is an abbreviation for y/l — f 2 . Tak- 
ing account of the self-consistency equation (11), we ob- 
tain the energy difference between the normal and super- 
conducting states: 

Fs ~ Fn (44) 



27rTA^(0) 



u»0 



f(f) , 9(9) 1 



2r+ 



One can check that this reduces to the known result for 
isotropic s-wave cases with or without pair breaking. 8 
This offers a straightforward way to calculate the the spe- 
cific heat near T c . The calculation, in general, is tedious 
because one has to keep track of terms up to A 4 oc St 2 . 
We consider only the case (A) = 0. 

Up to the forth order in A we have with the help of 
Eqs. (16) and (19): 

(A 2 ) 

v ; ■ 2 • (45) 



/ 



A 



2lo% 



+ 



-A s 

2r+w+ 
3A 4 A 2 (A 2 ) 



(46) 



2uj+ ' 8w 4 At + w\ ' 

where all u's are taken at T c . Substituting these in the 
energy difference we obtain: 

F s -F n 1 \ 

4 ^\ojI 2t+lo\ ) ' 



27rTAf(0) 



(47) 
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For large p+ one finds: 



i 



(3(fi 4 )-2)r 2 



3ttT 



(48) 



To complete the energy evaluation one needs *f>(T) 
which is obtained with the help of the self-consistency 
equation (27) and the expression (45) for /: 



^ 2 _47r 2 T c 2 (l-t) 



3(Q 4 > 



(49) 



Thus the energy difference between the normal and su- 
perconducting states reads: 



F s = 



8tt 4 N(0)i 



• k%T 2 c {T c - Tf 



(50) 



in common units. The specific heat jump at T c follows 



AC Ce Cjj 



16tt 4 /c 4 A^(0)r 2 
3ft 2 (3(ft 4 ) -2) 



(51) 



Within a weak coupling scheme, this result in a more 
general form has been obtained in Ref. 12. 

For the d-wave state on a cylindrical Fermi surface 
il = \72cos 2 20 and (fl 4 ) = 3/2 this gives: 



AC: 



32^ 4 4iV(0)r 2 ^ 



15ft 2 



(52) 



DISCUSSION 



Figure 1 is a compilation of data on the slopes H' c2 
for 1111 compounds with various dopants and, therefore, 
with various T^'s. An approximate scaling H' c2 oc T c is 
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FIG. 1: (Color online) The slopes of H c2 (T) near T c (the 
absolute values) for a few 1111 compounds. The data for the 
first three compounds in the legend are taken from Ref. 14; 
the remaining three points are taken from Ref. 15. The two 
right-most points are for H' c2ab of crystalline samples; the rest 
are for polycrystals, so that all points, in fact, reflect H' c2 ab . 



evident despite the fact that the compounds examined 
have T c 's varying from 6 to 46 K. From this data one 
estimates the slope of dH' c2 /dT c as « 0.2 T/K 2 . Then, 
the order of magnitude of the Fermi velocity follows from 
\dH' c2 /dT c \ ~ Tnfiokg/tfv 2 as v ~ 10 7 cm/s, a reasonable 
order that can be taken as yet another argument in favor 
of the picture presented. 

In Fig. 2 the data for the 122 family are collected. The 
same approximate scaling is seen. A considerable scatter 
of the data points might be caused by variety of rea- 
sons: different criteria in extracting H c2 from resistivity 
data, unavoidable uncertainties rooted in sample inhomo- 
geneities in determination of T c and the slopes of H c2 (T) 
near T c , possible differences in Fermi velocities and the 
order parameter anisotropies, to name a few. Moreover, 
the model employing only two scattering parameters for 
multi-band iron-pnictides is a far-reaching simplification, 
so that one can expect the model to work qualitatively at 
best. Nevertheless, the observed scaling seems remark- 
ably robust. One can take this as evidence in favor of a 
strong pair-breaking present in all compounds examined. 
It should be stressed again that for strongly anisotropic 
order parameters, (A) 0, the T c suppression (or the 
pair-breaking, which is the same) is caused by the com- 
bined effect of the transport and the spin-flip scattering. 

Figure 3 shows the specific heat jump measured in a 
number of compounds and reported in Ref. 1 . The "Ames 
scaling" AC oc suggested by Bud'ko, Ni, and Can- 
field is evident. Again, it is worth noting that only the 
combined rate p+ enters the coefficient in front of T c 3 of 
Eq. (51), so that the source of T c suppression is not nec- 
essarily the spin-flip AG pair-breaking. The ever present 
transport scattering suppresses T c as well, provided the 
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FIG. 2: (Color online) The slopes H c2 {T) near T c for a few 
122 iron-pnictides. The data are taken from: RbFe2As2 - 
Ref. 17, KFe 2 As 2 - Ref. 18, Ba .55Ko.45Fe 2 As2 - Ref. 16, the 
underdoped (ud) and overdoped (od) Ba(Fei_ x Co a ;)2Fe2As2 

- Ref. 3, Ba .6Ko.4Fe 2 As2 - Ref. 19, Sr .6Ko.4Fe 2 As2 - 
Ref. 20, Ba(Fe-Ni) 2 Fe 2 As2 and Ba(Fei_ 2 ,_ y Co a; Cuj / )2Fe2As2 

- Ref. 21. 
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FIG. 3: (Color online) The specific heat jump versus T c for a 
few 122 compounds shown on a log-log plot. The dashed line 
corresponds to AC oc T 3 . Most of the data are from Ref. 1; 
the new data points for mixed Co-Cu doping are shown by 
stars and taken by the same group, but have not been included 
in the original publication. 



order parameter is strongly anisotropic. This is presum- 
ably the case of iron-pnictides. 

One may wonder why the scaling H' c2 oc T c and AC oc 
seem to work across the whole class of iron pnictides 
for compounds with different couplings, F-surfaces etc. 
Clearly, the source of this scaling should be universal as- 
cross the phictides family of materials. The pair breaking 
is offered here as such an universal source. 

As for the apparent simplicity of the model used one 
should have in mind the often overlooked strength of the 
weak-coupling scheme: the model is formulated in terms 
of the measured critical temperature T c , in which the 
coupling constants and energy scales of the "glue bosons" 
are incorporated. 

Having succeeded in describing the "Ames scalings" 



just discussed, one can venture to a prediction: according 
to Eq. (41), tunneling experiments are likely to show the 
ratio of the apparent gap (the maximum position of the 
total density of states) to T c varying as 1/T C across the 
family of iron pnictides. 
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APPENDIX A: MATERIALS WITH (fi) / NEAR 

T c 

Interestingly enough, the behavior of the H c2 slopes as 
functions of T c turns out different if (fi) ^ 0. To see this, 
consider the expression for the coefficient A of Eq. (34). 
In terms of scattering times, it reads: 



.4 



2^T c r 2 



■kT c , 



2t_ 

T-L 



1 In- 



2r 4 



(Al) 



Since all r's are finite near the critical point where T c — > 
0, the term oc (f2) 2 is leading. Consider, e.g., a usual 
situation r << r m : 



.4 



(0) 2 ln2 
TrT r 



(A2) 



After simple algebra one obtains the slope of H c2c at T c : 

dH c2 , c = _Jo (») 2 ln4 

dT 2itt 2 T c (v 2 a W) + (ft) (v*Sl) ln(2r m /re 3 / 2 ) ' 

(A3) 

Thus, the slopes H' c2 oc 1/T C , the dependence opposite 
to that of the case (f2) = 0. 
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